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ABSTRACT 

We report the results of pulsar searches in ten globular clusters using the 
Robert C. Byrd Green Bank Telescope. One new binary millisecond pulsar has 
been discovered in NGC 5986 with P sp i n = 2.6 ms, P or b = 1.3 d, and a min- 
imum companion mass of 0.16 M . The companion is most likely a Helium 
white dwarf. Eight of the globular clusters we searched have central densities 
< 10 4 L Q pc -3 , making this a good sample for studying the pulsar content of 
low density clusters. We find no evidence for pulsars in clusters with very low 
densities (< 10 3 L Q pc -3 ), consistent with theoretical predictions. Null results 
in many of clusters we searched with moderate densities indicate that these sys- 
tems do not contain a bright MSP. Two clusters in particu lar, one with very low 



metallicity, stand in contrast to theoretical calculations bv llvanova et al.l (|20 08'l. 
We a lso find that three body exchange interaction rates calculated by (jphinney 



19961 ) seem to over predict the pulsar content in the clusters we studied. 



Subject headings: pulsars: individual (J1546— 3747A) — pulsars: general — globular 
clusters: individual (M92, NGC 288, NGC 2298, NGC 5897, NGC 5986, NGC 6981, 
NGC 7089, Pal 6, Pal 12, Terzan 9) — globular clusters: general 
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Introduction 



Over the past 23 years, 143 pulsars have been discovered in 27 globular clusters (GCs), 
the vast majority of which are true millisecond pulsars (MSPs) with P s < 10 njj]. Many 



of these MSPs are in exotic binary system s (e.g. iD'Amico et al.ll2001t iRansom et al.l 12005 
Hessels et al.ll2006t iFreire et al.1120071 120081 ). This stands in stark contrast to the population 
of pulsars found in the Galactic disk, which are dominated by isolated pulsars with longer 
spin periods and a small number of binary MSPs. These differences have been attributed to 
the high stellar densities in GC cores, where exchange interactions form mass-transferring 
binaries leading to the recycling of a "dead" pulsar. These exchange interactions are also 



responsible for the exotic binaries described above (see lCamilo fc Rasioll2005l ; lRansomll2008 
for reviews). The large number of exchange interactions in GCs can in turn be traced 
to the high stellar densities in cluster cores, and indeed, correlations between densi t y and 
the number of neutron stars in GCs have been known for some time. iPooley et al.l (120031 ) 
demonstrated that the core interaction rate, T c oc pl' 5 r 2 , is a good predictor of the number 
of low-mass X-ray binaries, which are thought to be the progenitor systems to many MSPs. 
The value of T c /D 2 , where D is the distance and accounts for flux losses, has been used 
to pick out the mos t promising clusters for pulsar searche s with great success. Terzan 5 
j Ransom et al.l 120051 ; iHessels et al.l 120061 ) and 47 Tucanae ( iCamilo et al.l |2000e IFreire et al. 
20031 ) are the two richest clusters for MSPs and have the highest value of T c /D 2 . Recent 
simulati ons also show a str ong dependence on core density for the number of observable 
pulsars (llvanova et al.l 120081 ) . 



Nonetheless, the re are exceptions to this rule that have been known for some time. 
Kulkarni et al.l (ll99ll ) discovered two pulsars in the low density GCs M13 and M53. Sub- 
sequent searches uncovered a substantial population of MSPs in GCs with p c ~ 10 3 - 



10 4 L pc 3 . The to tals now stand at five MSPs in M5 (IWolszczan et al. 



2007), five in M13 flKulkarni et al.l Il991 



(IHessels et al.ll2007f ). t wo in NGC 6749 (IHessels et al 



Anderson 11993: 



1989; Hessels et al. 



200 



Hessels et all 120071). four in M3 
3), one in M53 ( Kulkarni et al. 



19911 ). and one in M71 (IHessels et al.ll2007l ). These MSPs account for about 13% of all known 
GC pulsars, despite the clusters contributing only 3 % to the total T c of a ll GCs with pulsars. 
However, most of these pulsars were discovered by IHessels et al.l (120071 ) using Arecibo, and 
without the exquisite sensitivity of this tele scope, only a handful of p ulsars would haye been 
found. Shortly after iKulkarni et al.l (Il99ll ) announced their results, iJohnston et al.l (119921 ) 
argued that pulsar formation in low density systems was enhanced relative to high density 
systems, finding N psy: oc p^ 5 (though this came prior to many fruitful searches in higher 



^ee http://www.naic.edu/~pfreire/GCpsr.html for an up-to-date list 
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density clusters). ISigurdsson h Phinneyl ( 119951 ) suggested that three-body encounters with 
primordial binaries could give rise to a small population of pulsars even in low-density GCs. 

We have searched ten clusters using the Robert C. Byrd Green Bank Telescope (GBT). 
We selected GCs that had not been searched using the GBT, focusing primarily on clusters 
that had low central densities and large Galactic latitudes (to avoid oversubscribed LST 
ranges). Five clusters (NGC 288, NGC 2 298, NGC 6981 NG C 7089, and Pal 12) were 
searched using the Parkes radio telescope (IPossenti et al.l 120051 ) . but we ensured that our 
searches were at least twice as sensitive. Since eight of the clusters have p c < 10 4 L pc~ 3 , 
this a good sample of deeply searched low density systems. We have discovered one new 
binary MSP in NGC 5986, and we discuss the implications of non-detections in the other 
clusters. 



2. Survey Details 

All our searches were car ried out using the G BT and the Green Bank Ultimate Pulsar 



Processor (GUPPI) back-end (jDuPlain et al.ll2008l ). Observations were made at a central fre- 



quency of 2 GHz, with 800 MHz of bandwidth broken into 2048 frequency channels, although 
persistent radio frequency interference (RFI) reduces the usable bandwidth to 600 MHz. We 
used sampling times of 40.96-64 /xs. Total system temperatures were typically 24-28 K. 
The contribut i on fro m the Galactic background was estimated by scaling the values from 



Haslam et al.l (119821 ) to 2 GHz assuming a spectral index of —2.6, though this was usually 
< 1 K except for clusters at very low Galactic latitudes. The limiting flux density of each 
search was taken to be 



ge^tot / w 

Gy/n pol Aut int V W-P 

where (3 = 1.3 accounts for digitization losses and £ = 10 is our limiting S/N. T tot is the 
total sky and system temperature, G = 1.9 is the telescope gain, n po \ = 2 is the number of 
summed polarizations, Au is the usable bandwidth, and t- mt is the integration time. P is the 
pulse period and W is the observed pulse width, which is a combination of the intrinsic pulse 
width, instrumental time resolution, and broadening of the pulse due to dispersive smearing 
and scattering. This calculation does not take into account sensitivity losses due to RFI or 
Doppler smearing of the pulse due to binary acceleration. However, we observed very little 
RFI in our data, and since we performed acceleration searches (see below), only pulsars in 
very tight binaries should have dropped below our detection threshold. Integration times 
and approximate limiting flux densities can be found in Table (U along with some other 
properties of each cluster. 
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To test the accuracy of our limiting flux density estimates, we calculated the integration 
time necessary to just detect the lone pulsar we discovered in NGC 5986, and then blindly 
searched this portion of our original data. The pulsar was indeed detected with with an 
appropriate S/N. We are therefore confident that our reported limits are accurate. 

Data were reduced using the PRESTcl software suite (IRansom et al.ll2002l ). After remov- 
ing RFI, de-dispersed time series were created, starting with a minimum dispersion measure 
(DM) of 10 pc cm" 3 and increasing in st eps of 0.1 pccm -3 to twice the DM predicted by 
the NE2001 model (ICordes fc Lazid 120021 ) . Acceleration techniques were used to search for 
periodic signals from isolated and binary MSPs and single pulse searches were used to look 
for bright transient emission. To search for highly accelerated pulsars, we performed ac- 
celeration searches on ~ 10-20 minute subsets of each observation. We searched up to a 
maximum acceleration of 800 Fourier bins in all acceleration searches. Candidate pulsars 
were visually inspected and grouped into likely pulsars or random noise and RFI. We were 
able to observe some clusters twice, which allowed us to confirm or reject candidate pulsars 
quickly. 



3. The Binary Pulsar in NGC 5986 

We discovered one new pulsar in NGC 5986, J1546-3747A, hereafter NGC 5986A (Fig- 
ure [TJ. The pulsar was discovered in our acceleration searches with an acceleration of —26 
Fourier bins, or —1.1 ms -2 , indicating that it is in a binary system. We conducted three 
follow-up observations to determine the orbit of NGC 5986A and to search for more pulsars 
using GUPPI in a coherent de-dispersion mode at 820 MHz, 1.5 GHz, and 2 GHz. The pulsar 
was detected in all three follow-up observations. We searched the 820 MHz and 2 GHz ob- 
servations, but found no new pulsars. The 1.5 GHz observations were heavily contaminated 
with RFI, making a blind search impractical. 

While these observations are not sufficient to obtain a full timing solution for NGC 5986A, 
they have allowed us to accurately determine the orbital parameters (see Tableland Fig. [2J). 
The orbit is circular, with a period of ~ 1.3 days and projected semi-major axis of ~ 0.6 R . 
The companion mass limit is 0.16 M assuming a pulsar mass of 1.4 M Q . Given a distribu- 
tion of orbital inclinations that is flat in cost, an edge-on orbit is statistically most likely, 
so the true companion mass is probably close to this limit. These binary characteristics 
are similar to those of MSPs with Helium white dwarf companions that are found in both 
GCs and the Galactic disk. A main-sequence companion seems unlikely, as such systems 



http : //www . cv . nrao . edu/~sransom/presto/ 
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have only been found in dense GCs ([Edmonds et al.ll2002l ; iPossenti et al.ll2003l ; iHessels et al. 
20061 ). Furthermore, we see no evidence for eclipses, although we lack coverage near conjunc- 
tion when eclipses would be most likely. Systems similar to NGC 5986A have been found 
in other GCs of comparable densityE). If a precise position can be measured through future 
pulsar timing, a search for an optical counterpart may be feasible. 

Rough flux density estimates were obtained by assuming that the off-pulse RMS noise 
levels are described by the radiometer equation, 



a 



tot 



Gy/n pol Avt 



(2) 



mt 



We find S u = 21, 27, and 45 /iJy at 2 GHz, 1.5 GHz, and 820 MHz, respectively, which 
yields a spectral index of a ~ —0.9, flatter than the canonical value of —1.7. 



4. The Pulsar Content of Low-Density Clusters 



We searched ten clusters for pulsars, five of which have 10 3 < p f <4x 10 4 pc 3 , but 



found only one MSP. This stands in contrast to previous surveys (iHessels et al.l 120071 . and 
references therein, hereafter IHRS07I ) that found several pulsars in low density clusters. We 
therefore consider what factors may have led to our null results. 

Fig. Oshows the approximate limiting luminosity of our searches, where L min = S min D 2 , 
along with the 1.4 GHz pseudo-luminosities of pulsars with reliable flux density measure- 
ments from IHRS07I . scaled to 2 GHz u sing a s pectral index of —1.7 (NGC 5986A is also 
shown). It is immediately obvious that IHRS07I were more sensitive than o ur searches. We 
can attribute this to two factors. The first is that all the GCs from IHRS07I had D < 8 kpc, 
with the exception of M53 (D = 17.8 kpc). By comparison, only two clus ters in o ur survey 
have D < 8 kpc, while most have D > 10 kpc. The second factor is that IHRS07I were able 
to reach lower limiting flux densities using Arecibo than we could with the GBT. 

Despite these factors, 3/10 of ou r search es should have been sensitive to the brightest 
~ 53% of a population similar to the IHRS07I pulsars, and 8/10 of our searches should have 
been sensitive to the brightest two pulsars. Furthermore, four out of five of the clusters 
with p c > 10 3 L pc -3 were among the five most sensitively searched GCs in our sample. In 
other words, we achieved good sensitivity in the most promising clusters. As discussed in 
£J2J we are confident that our limiting flux densities are accurate unless there are pulsars in 



"http : //www . naic . edu/~pf reire/ GCpsr . html 
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extremely tight binaries. However, since IHRS07I used the same search procedure as we have, 
they suffered from the same bias, so we do not believe this is a good exp lanation for the lack 
of pulsars in our sample. While it is possible that some of the IHRS07I pulsars could have a 
steeper spectral index than the assumed value of —1.7, we would also expect some to have 
flatter spectra, making them easier to detect at 2 GHz. It therefore seems likely that most 
of the clusters we searched lack a bright MSP. 



We now turn to the theoretical results of llvanova et al.l ( 120081 ). who predict the number 
of neutron stars in different types of GCs through Monte Carlo simulations. They simulated 
a variety of clusters, including those with lo g (/? r /pc~ 3 ) = 3 , 4,5,6 . We began by empirically 



fitting a simple exponential to the results of llvanova et al.l (120081 ) of the form 

3 logp c _|_ £ 



n 



psrs 



a e 



(3) 



with a = 0. 041 and b = —0.82, a nd where n psrs is the number of pulsars formed per 2 x 10 5 M Q 
(as per the llvanova et al.l (120081 ) simulations). We explored other functional forms but chose 
this one for its goodness of fit and simplicity. We then calculated the expected number of 
pulsars in our clusters based upon their central densities and masses^. Next, we calculated 
how many pulsars would lie above our limiting flux densities (at 2 ms), given a power-law 
distribution, dN(L) oc L" 1 dL. We assumed a maximum luminosity equal to the brightest 
cluster pulsars, (about 250 mJykpc 2 at 2 GHz) and a minimum luminosity of 0.16 mJykpc 2 
(obtained by scaling the typically assumed lower limit of 0.3 mJykpc 2 at 1.4 GHz to 2 GHz). 
We calculated upper and lower estimates in the same way, fitting to the maxim um and 
minimum expected pulsars as defined by the errors quoted in llvanova et al.l (120081 ) . 



Based on these calculations, we expect that about nine pulsars should have been detected 
in our sample, with upper and lower estimates of 16 and three pulsars, respectively. Our 
results are not very sensitive to our choices of L max and L min since our limiting luminosities are 
above the L min cutoff. Choosing a flatter slope for the luminosity distribution only increases 
the number of potentially observable pulsars. As a chec k on our method, we performed the 
same analysis using all the clusters searched by IHRS07I . M15 (a dense, massive cluster with 
eight known pulsars) was an outlier in these calculations, with far more pulsars predicted 
than are obser ved; when it is excluded from the analysis, we predict 6^g 5 observable pulsars 
in the IHRS07I sample. Twenty have been discovered (excluding the eight in M15), which is 
consistent with this upper limit. 



Phinneyl (119961 ) also gave birth rates for pulsars in GCs for different formation mecha- 



nisms, the most relevant in our density regime being three-body exchanges. The birth rates 



4 Masses were taken from photometric models by O. Y. Gnedin; 
http : //www. astro . lsa.umich.edu/~ogriedin/gc/vesc . dat 
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are given per neutron star, so we once again use the results of llvanova et al.l (120081 ) for an 
estimate of the total neutron star content of our clusters. We also assume an MSP lifetime 
of ~ 10 Gyr. Based on these theoretical models, we should have been sensitive to ~ 1-2 
pulsars each in NGC2298, Terzan 9, and NGC6981, about ten in Pal 6, about two dozen in 
NGC5986, ~ 50 in M92, and ~ 70 in NGC7089. 



Discussion 



We p redict that 9^ pulsa rs should have been observable in our sample based on the 
results of llvanova et al.l ( 120081 ). This prediction seems high given our single detection in 
NGC 5986 (we note that this cluster was pred icted to contain one observable pulsar based 
on the high estimates of llvanova et al.l (120081 ) ) . In addition, nearly all of the potentially 
observable MSPs come from M92 and NGC 7089. This assertion is supported by the fact that 
these two clusters have the highest value of T c / D 2 in our sample. As can be seen from Fig. 
El M92 in particular was searched fairly deeply, so this null result is particularly interesting. 
One parameter that was not taken into account in our simulations was metallicity. Surveys 
of Galactic and extra-galactic GCs h ave found ~ 3 times more bright LMXB s in metal-rich 



2004 




Kundu et al. 


2002 





clust e rs (defined by [Fe/Hl < -1) (ISivakoff et al.l 120071 : IJordan et al.l 120041 : iMinniti et al. 



effect may also manifest itself in the MSP population. We note that eight of the clusters 
we studied have [Fe/H] < — 1, and that M92 in particular has the second lowest metallicity 
of any Milky Way GC — the most metal poor is M15, which was greatly under-populated in 
MSPs based on the calculations discussed in S]H While a full treatment of metallicity effects 
is beyond the scope of this paper, our results may suggest that MSPs are less likely to be 
found in low metallicity clusters. 



We find that the results of iPhinneyl (119961 ) drastically over-predict the number of MSPs 
that should be observable in our study. One possible explanation for this discrepancy is 
that our assumed MSP lifetime of 10 Gyr is too long. An order of magnitude decrease in 
the lifetime would yield numbers that are closer to our observed results. However, based 
on observed spin down rates, MSPs should be very long lived. The other par ameters in our 



calcul ations were the number of primordial neutron stars (which we based off of llvanova et al. 



2008m) and the rate of MSP formation through exchange encounters calculated by iPhinney 
( 119961 ) . It therefor e seems t h at on e or both of these inputs is too large for clusters of low to 



moderate density. IPhinneyl (119961 ) over predicts the number of observable MSPs in higher 
density clusters as well, so it seems that these rates may not be applicable when trying to 
predict the number of MSPs in a given cluster. 
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We thus draw the following conclusions. GC pulsar searches are still sensitivity limited, 
and this necessarily makes any statements about the total population of cluster pulsars a 
matter of extrapolation. Nevertheless, our results show no evidence for a population of 
pulsars in v ery low density syst ems (< 10 3 L pc -3 ), which is consistent with theoretical 



predictions (jlvanova et al.l 120081 ) . Even given our sensitivity limits, though, it appears that 



cl usters w ith densities 10 3 -10 4 L pc 3 are not as efficient at forming MSPs as the results 



of IHRS07I imply. A full treatment of the effects of metallicity may shed more light on this 
discrepancy. We find evidence that low (and high) density GCs are either not as efficient at 
forming neutro n stars or recyc ling; them into MSPs through exchange interactions as earlier 



calculations by iPhinneyl (jl996[ ) imply. 
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Table 1. Globular Clusters Included in Survey 





t 

(deg) 


Z. 


(deg) 


u 

(kpc) 


-rreuicteu JJM 
(pc cm -3 ) 


r c 

(arcmin) 


log Pc 

(L© pc ) 


t / n2 

to? a, \ 

(7o ) 


[le/HJ 


^obs 


Approximate S v min 

(/~Jy) 


NC!C! 988 




—89 38 


8 9 


98 


1 35 


1 78 

1. IO 


n D3 


— 1 39 


n 95 


1 1 


NGC 2298 


245.63 


-16.00 


10.8 


85 


0.31 


2.90 


0.08 


-1.92 


3.1 


6 


NGC 5897 


342.95 


+30.29 


12.5 


61 


1.40 


1.53 


0.01 


-1.90 


1.3 


10 


NGC 5986 


337.02 


+13.27 


10.4 


92.2 


0.47 


3.41 


1.1 


-1.59 


2.2 


7 


M92 


68.34 


+34.86 


8.3 


43 


0.26 


4.30 


7.0 


-2.31 


1.4 


9 


Pal 6 


2.10 


1.78 


5.8 


397 


0.66 


3.46 


2.5 


-0.91 


1.2 


15 


Terzan 9 


3.61 


-1.99 


7.1 


422 


0.03 


4.42 


0.14 


-1.05 


2.5 


10 


NGC 6981 


35.16 


-32.68 


17.0 


55 


0.46 


2.38 


0.03 


-1.42 


1.0 


11 


NGC 7089 


53.37 


-35.77 


11.5 


46 


0.32 


4.00 


3.8 


-1.65 


1.0 


11 


Pal 12 


30.51 


-47.68 


19.0 


40 


0.02 


3.64 


0.004 


-0.85 


1.0 


11 


47 Tucanae b 


305.89 


-44.89 


4.5 


24.4 


0.36 


4.88 


100 


-0.72 







a Normalized to 47 Tucanae. 

b 47 Tucanae has been included so that the properties of our sample can be compared to a rich GC. 



Note. Cluster parameters are taken from the 2010 version of the Harris catalog 

(http://physwww.physics.mcmaster.ca/~harris/mwgc.dat). DM estimates for all clusters except NGC 5986 and 47 
Tucanae are from the NE2001 model (ICordes Lazioll2002l ). 



- 12 - 



Table 2. Timing Derived Parameters of NGC 5986 



Parameter 



Value 



P (ms) 

Reference Epoch (MJD) 
DM (pc cm" 3 ) 

asmi/ c (s) 
T (MJD) 

e 

u (deg) 

Mass Function (M ) 



-^c,min 
NtqAs 



(M G 



0; 



2.6056722466(1) 
51900 
92.17(4) 
1.3467116(2) 
1.38325(4) 
55355.67467(2) 
< 8 x 10~ 6 

0.0016 
0.16 
79 



"Assuming M psr = 1.4 M 



Note. All parameters assume that 

the pulsar is at the cluste r center, a = 
15:46:03.44, 5 = -37:47:10.1 Jshawl & White 



1986J). The eccentricity limit was calculated 
according to en m = St (asmi/ c) -1 w here St is 
our timing precision ( lPhinneylll992l ). The or- 
bital solution was obtained using the DE405 
Solar System ephemeris and the UTC(NIST) 
time standard. 
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Pulse Phase 

Fig. 1. — The flux-calibrated pulse profile of NGC 5986A at an observing frequency of 2 GHz 
(S v = 21 fjjy). 



- 14 - 



Days 




Orbital Phase 



Fig. 2. — The Doppler modulated period of NGC 5986A, clearly showing the effects of 
binary acceleration and a circular orbit. The thick lines show the observed period during 
our observations, and the thin line is the predicted period based on our orbital solution. 
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Fig. 3. — Comparison of the limiting luminosity of our s earches and the luminosities of 
pulsars with reliable flux densities from lHessels et al.l (120071 ). as well as NGC 5986A. Lumi- 
nosities of the known MSPs were scaled from 1.4 GHz to 2 GHz using a spectral index of —1.7. 
Solid lines indicate clusters with p c > 10 4 L pc -3 , dashed lines indicate 10 3 -10 4 L pc~ 3 , 
and dot-dashed lines indicate < 10 3 L pc~ 3 . 



